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Abstract. We first show that a continuous function / is nonnegative on a 
closed set K C R n if and only if (countably many) moment matrices of some 
signed measure du = fd/i with supp fj, = K, are all positive semidefinite (if 
K is compact fi is an arbitrary finite Borel measure with supp fi = K). In 
particular, we obtain a convergent explicit hierarchy of semidefinite (outer) 
approximations with no lifting, of the cone of nonnegative polynomials of de- 
gree at most d. Wen used in polynomial optimization on certain simple closed 
sets K (like e.g., the whole space R n , the positive orthant, a box, a simplex, 
or the vertices of the hypercube), it provides a nonincreasing sequence of up- 
per bounds which converges to the global minimum by solving a hierarchy of 
semidefinite programs with only one variable (in fact, a generalized eigenvalue 
problem). In the compact case, this convergent sequence of upper bounds 
complements the convergent sequence of lower bounds obtained by solving a 
hierarchy of semidefinite relaxations as in e.g. |12| . 



1. Introduction 

This paper is concerned with a concrete characterization of continuous functions 
that are nonnegative on a closed set K C R n and its application for optimization 
purposes. By concrete we mean a systematic procedure, e.g. a numerical test that 
can be implemented by an algorithm, at least in some interesting cases. For polyno- 
mials, Stengle's Nichtnegativstellensatz [22] provides a certificate of nonnegativity 
(or absence of nonnegativity) on a semi-algebraic set. Moreover, in principle, this 
certificate can be obtained by solving a single semidefinite program (although the 
size of this semidefinite program is far beyond the capabilities of today's comput- 
ers). Similarly, for compact basic semi-algebraic sets, Schmudgen's and Putinar's 
Positivstellensatze [SU] [18] provide certificates of strict positivity that can be ob- 
tained by solving finitely many semidefinite programs (of increasing size) . Exten- 
sions of those certificates to some algebras of non-polynomial functions have been 
recently proposed in Lasserre and Putinar [14] and in Marshall and Netzer [16] . 
However, and to the best of our knowledge, there is still no hierarchy of explicit 
(outer or inner) semidefinite approximations (with or without lifting) of the cone 
of polynomials nonnegative on a closed set K, except if K is compact and basic 
semi-algebraic (in which case outer approximations exist). Another exception is 
the convex cone of quadratic forms nonnegative on K = K" for which inner and 
outer approximations are available; see e.g. Anstreicher and Burer pQ, and Diir [7]. 
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Contribution: In this paper, we present a different approach based on a new 
(at least to the best of our knowledge) and simple characterization of continuous 
functions that are nonnegative on a closed set K C M™. This characterization 
involves a single (but known) measure fi with supp /i = K, and sums of squares of 
polynomials. Namely, our contribution is twofold: 

(a) We first show that a continuous function / is nonnegative on a closed set 
K C R n if and only if J h 2 fdfi is nonnegative for all polynomials h g K[x], where 
fi is a finite Borel measurqj with supp /x = K. The measure n is arbitrary if K is 
compact. If K is not compact then one may choose for /i the finite Borel measure: 

- d/i = exp(— \xi\)d(f if / is a polynomial, and 

- dfJ, = (1 + f 2 )^ 1 exp(— J2 t \xi\)d<p, if / is not a polynomial, 

where ip is any finite Borel measure with support exactly K. But many other 
choices are possible. 

Equivalently, / is nonnegative on K if and only if every element of the countable 
family T of moment matrices associated with the signed Borel measure /d/x, is pos- 
itive semidefinite. The absence of nonnegativity on K can be certified by exhibiting 
a polynomial h £ R[x] such that J h 2 fd[i < 0, or equivalently, when some moment 
matrix in the family T is not positive semidefinite. And so, interestingly, as for 
nonnegativity or positivity, our certificate for absence of nonnegativity is also in 
terms of sums of squares. When / is a polynomial, these moment matrices are eas- 
ily obtained from the moments of /x and this criterion for absence of nonnegativity 
complements Stengle's Nichtnegativstellensatz [35] (which provides a certificate of 
nonnegativity on a semi-algebraic set K) or Schmiidgen and Putinar's Positivstel- 
lensatze 20, 18] (for certificates of strict positivity on compact basic semi-algebraic 
sets). At last but not least, we obtain a convergent explicit hierarchy of semidefinite 
(outer) approximations with no lifting, of the cone Cd of nonnegative polynomials 
of degree at most 2d. That is, we obtain a nested sequence C° D ■ ■ - C^ D ■ ■ ■ D Cd 
such that each C\ is a spectrahedron defined solely in terms of the vector of coeffi- 
cients of the polynomial, with no additional variable (i.e., no projection is needed). 
Similar explicit hierarchies can be obtained for the cone of polynomials nonnegative 
on a closed set K (neither necessarily basic semi- algebraic nor compact), provided 
that all moments of an appropriate measure /i (with support exactly K) can be 
obtained. To the best of our knowledge, this is first result of this kind. 

(b) As a potential application, we consider the problem of computing the global 
minimum /* of a polynomial / on a closed set K, a notoriously difficult problem. 
In nonlinear programming, a sequence of upper bounds on /* is usually obtained 
from a sequence of feasible points (x&) C K, e.g., via some (local) minimization 
algorithm. But it is important to emphasize that for non convex problems, provid- 
ing a sequence of upper bounds (/(xfc)), k £ N, that converges to /* is in general 
impossible, unless one computes points on a grid whose mesh size tends to zero. 

We consider the case where K C l n is a closed set for which one may compute 
all moments of a measure /x with supp /x = K. Typical examples of such sets are 
e.g. K = M™ or K = M™ in the non compact case and a box, a ball, an ellipsoid, a 
simplex, or the vertices of an hypercube (or hyper rectangle) in the compact case. 



A finite Borel measure /i on R n is a nonnegative set function defined on the Borel tr-algebra 
of M n (i.e., the cr-algebra generated by the open sets), such that /x(0) = 0, ^t(R n ) < oo, and 
^(USi Ei) = J^Si M^i) f° r an y collection of disjoint measurable sets Bj. Its support (denoted 
supp fi) is the smallest closed set K such that \ K) = 0; see e.g. Royden [19] . 
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We then provide a hierarchy of semidefinite programs (with only one variable!) 
whose optimal values form a monotone nonincreasing sequence of upper bounds 
which converges to the global minimum /*. In fact, each semidefinite program is 
a very specific one as it reduces to solving a generalized eigenvalue problem for a 
pair of real symmetric matrices.. (Therefore, for efficiency one may use specialized 
software packages instead of a SDP solver.) However, the convergence to /* is in 
general only asymptotic and not finite (except when K is a discrete set in which case 
finite convergence takes place) . This is in contrast with the hierarchy of semidefinite 
relaxations defined in Lasserre |12[ 113] which provide a nondecreasing sequence of 
lower bounds that also converges to /*, and very often in finitely many steps. 
Hence, for compact basic semi-algebraic sets these two convergent hierarchies of 
upper and lower bounds complement each other and permit to locate the global 
minimum /* in smaller and smaller intervals. 

Notice that convergence of the hierarchy of convex relaxations in [12] is guar- 
anteed only for compact basic semi-algebraic sets, whereas for the new hierarchy 
of upper bounds, the only requirement on K is to know all moments of a measure 
fj, with supp/Lt = K. On the other hand, in general computing such moments is 
possible only for relatively simple (but not necessarily compact nor semi-algebraic) 
sets. 

At last but not least, the nonincreasing sequence of upper bounds converges to /* 
even if /* is not attained, which when K = R™, could provide an alternative and/or 
a complement to the hierarchy of convex relaxations provided in Schweighofer |21j 
(based on gradient tentacles) and in Ha and Pham [8] (based on the truncated 
tangency variety), which both provide again a monotone sequence of lower bounds. 

Finally, we also give a very simple interpretation of the hierarchy of dual semi- 
definite programs, which provides some information on the location of global min- 
imizers. 



A Borel measure on M n is understood as a positive Borel measure, i.e., a non- 
negative set function /j, on the Borel tr-algebra B (i.e., the u-algebra generated by 
the open sets of R") such that /u(0) = 0, and with the countably additive property 



for any collection of disjoint measurable sets C B; see e.g. Royden [19l pp. 
253-254]. 

Let R[x] be the ring of polynomials in the variables x = (xi, . . . , x n ), and S[x] C 
R[x] its subset of polynomials that are sums of squares (s.o.s.). Denote by R[x]d C 
R[x] the vector space of polynomials of degree at most d, which forms a vector space 
of dimension s(d) = with e.g., the usual canonical basis (x Q ) of monomials. 

Also, denote by E[x]d C S[x] the convex cone of s.o.s. polynomials of degree at 
most 2d. If / € R[x]d, write /(x) = X^aeW /a xQ m the canonical basis and denote 
by f = (f a ) 6 R s ( rf ) its vector of coefficients. Let S n denotes the vector space of 
p x p real symmetric matrices. For a matrix A £ S p the notation A y (resp. 
A >- 0) stands for A is positive semidefinite (resp. definite). 



2. Notation, definitions and preliminary results 
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Moment matrix. With y = (y a ) being a sequence indexed in the canonical basis 
(x Q ) of R[x], let L y : R[x] -» M be the linear functional 

/ (=£/ a x a ) ^ L y (f) = £/ a y a , 

a a 

and let Md(y) be the symmetric matrix with rows and columns indexed in the 
canonical basis (x Q ), and defined by: 

(2.1) M d (y)(a,f3) := L y (x«+^) = y a+0 , a, e N% 

with := {a e N" : |a| (= «») < 4- 

If y has a representing measure fj,, i.e., if y a = Jx tt d/i for every a € N™, then 

(f,M d (y)f) = y/(x) 2 ^(x) > 0, V/GR[x] d , 

and so Md(y) ^0. A measure /x is said to be moment determinate if there is no 
other measure with same moments. In particular, and as an easy consequence of the 
Stone- Weierstrass theorem, every measure with compact support is determinat^. 

Not every sequence y satisfying Md(y) >: 0, d € N, has a representing measure. 
However: 

Proposition 2.1 (Berg [3]). Let y — (y a ) be such that Md(y) el 0, for every 
deN. Then: 

(a) The sequence y has a representing measure whose support is contained in the 
ball [—a,a] n if there exists a,c > such that \y a \ < ca) a > for every a G N™. 

(b) The sequence y has a representing measure \x on R" if 

oo 

(2.2) £ L y (xT)- 1/2t = +oo, Vi = 1, . . • , n. 

In addition, in both cases (a) and (b) the measure \i is moment determinate. 

Condition (b) is an extension to the multivariate case of Carleman's condition 
in the univariate case and is due to Nussbaum [17]. For more details see e.g. Berg 
[3] and/or Maserick and Berg [11] . 

Localizing matrix. Similarly, with y = (y a ) and / G R[x] written 

xi->/w = ^ x7 ' 

let M<j(/y) be the symmetric matrix with rows and columns indexed in the canon- 
ical basis (x Q ), and defined by: 

(2.3) M d (/y)(a,/3) := L y (/(x)x«+^) = ]T / 7 2/ Q+/3+7 , Va,/3eNS- 

7 

If y has a representing measure fj., then (g, M^/ y)g) = J g 2 fd[i, and so if /x is 
supported on the set {x : /(x) > 0}, then M^(/ y) y for all d = 0, 1, . . . because 

(2.4) (g,M d (/y)g) - /" s(x) 2 /(x)d/i(x) > 0, Vg e K[x] d . 

2 To see this note that (a) two measures /^i , /i2 on a compact set K C R" are identical if and 
only if J K /d/ii = J K /rf/i2 for all continuous functions / on K, and (b) by Stone- Weierstrass, the 
polynomials are dense in the space of continuous functions for the sup-norm. 
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3. NONNEGATIVITY ON CLOSED SETS 

Recall that if X is a separable metric space with Borel <7-field B, the support 
supp /j, of a Borel measure /i on X is the (unique) smallest closed set B £ B such 
that n{X \ B) = 0. Given a Borel measure fi on R™ and a measurable function 
/ : 1™ — > R, the mapping B i-> v(B) := J B fdfx, B £ B, defines a set function on B. 
If / is nonnegative then v is a Borel measure (which is finite if / is /i-integrable); 
see e.g. Royden [19j p. 276 and p. 408]. If / is not nonnegative then setting 
B\ := {x : /(x) > 0} and B 2 := {x : /(x) < 0}, the set function v can be written 
as the difference 

(3.1) v{B) = Vl (B)-u 2 (B), B e B, 
of the two positive Borel measures v\ , v 2 defined by 

(3.2) vi{B)=\ fdn, v 2 (B) = - [ fdn, VBeB. 

J Bins J B 2 nB 

Then v is a signed Borel measure provided that either Vx(Bx) or ^2(^2) is finite; 
see e.g. Royden [THl p. 271]. We first provide the following auxiliary result which 
is also of self-interest. 

Lemma 3.1. Let X be a separable metric space, K C X a closed set, and /i a Borel 
measure on X with supp/i = K. A continuous function f : X — ► R is nonnegative 
on K if and only if the set function B > v{B) = J KnB fdfj,, B 6 B, is a positive 
measure. 

Proof. The only if part is straightforward. For the if part, if v is a positive measure 
then /(x) > for ^(-almost all x € K. That is, there is a Borel set G C K such 
that fi(G) = and /(x) > on K \ G. Indeed, otherwise suppose that there 
exists a Borel set Bq with fi(Bo) > and / < on Bo; then one would get the 
contradiction that v is not positive because v(Bq) = j Bg fd/j, < 0. In fact, / is called 
the Radon-Nikodym derivative of v with respect to fi; see Royden pjH Theorem 23, 

p. 276]. 

Next, observe that K \ G C K and fi(K \ G) = /x(K). Therefore K \ G = K 
because supp/i(= K) is the unique smallest closed set such that /i(X \ K) = 0. 
Hence, let x G K be fixed, arbitrary. As K = K \ G, there is a sequence (x/.) C 
K \ G, k G N, with x,t — > x as k —> 00. But since / is continuous and /(x^) > 
for every k € N, we obtain the desired result /(x) > 0. □ 

Lemma lBTTI itself (of which we have not been able to find a trace in the literature) 
is a characterization of nonnegativity on K for a continuous function / on X. 
However, one goal of this paper is to provide a more concrete characterization. To 
do so we first consider the case of a compact set K C R™. 

3.1. The compact case. Let K be a compact subset of K". For simplicity, and 
with no loss of generality, we may and will assume that K C [— 1, l] n . 

Theorem 3.2. Let K C [—1, 1]™ be compact and let /j, be an arbitrary, fixed, finite 
Borel measure on K with supp \x = K, and with vector of moment y = (y a ), ot S N n . 
Let f be a continuous function on R™. Then: 
(a) / is nonnegative on K if and only if 

(3.3) f g 2 fdfi>0, VjeR[4 
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or, equivalently, if and only if 

(3.4) M d (z) >r 0, d=0,l,... 

where z = (z a ), a € N™ , wii/i z Q = Jx°/(x)rf/i(x), and im£/i M^(z) as in R2.1]) . 

If in addition f £ K[x] then {3.4\) reads M,j(/y) h 0, d — 0,1,..., where 
Mrf(/y) ?s i/ie localizing matrix defined in H2.3\) . 

(b) // m addition to be continuous, f is also concave on K, then f is nonnegative 
on the convex hull co (K) of K if and only if i3.3\) holds. 

Proof. The only if part is straightforward. Indeed, if / > on K then K C {x : 
/(x) > 0} and so for any finite Borel measure /ionK, J K g 2 fd[i > for every 
g £ R[x]. Next, if / is concave and / > on co (K) then / > on K and so the 
"only if part of (b) also follows. 

If part. The set function v(B) = J B fdfi, B £ B, can be written as the difference 
v = V2 — V2 of the two positive finite Borel measures v\, described in (|3.ip - (|3.2l) . 
where Bi :={xeK: /(x) > 0} and B 2 :={xeK: /(x) < 0}. As K is compact 
and / is continuous, both v\, v% are finite, and so v is a finite signed Borel measure; 
see Royden [19j p. 271]. In view of Lemma [3.11 it suffices to prove that in fact v 
is a finite and positive Borel measure. So let z = (z a ), a £ N", be the sequence 
defined by: 

(3.5) z a = f x Q aV(x) := f x a f{x)d f j,(x), Va £ N". 

Every z a , a £ N", is finite because K is compact and / is continuous. So the 
condition 

/ 5 (x) 2 /(x)dM(x)>0, V/£% 

reads (g,M d (z)g) > for all g £ W (d \ that is, M d (z) > 0, where M d (z) is the 
moment matrix defined in (|2.1|) . And so (|3.3|) implies Md(z) > for every del. 
Moreover, as K C [—1,1]™, 

\z a \ < c := f \f\dfi, VaeN". 

Hence, by Proposition 12.11 z is the moment sequence of a finite (positive) Borel 
measure i[> on [—1, 1]™, that is, as supp^ C K C [—1, 1]", 

(3.6) / x Q aV(x) = / x a #(x), VaeN". 
J [-1,1]" J[-lA] n 

But then using (|3.ip and (|3.6p yields 

/ x Q aVi(x) = / x a d{v 2 + V)(x), VaeN", 

which in turn implies v\ = v 2 + tp because measures on a compact set are determi- 
nate. Next, this implies ip — v\ — V2 (— v) and so v is a positive Borel measure on 
K. Hence by Lemma T3.ll /(x) > on K. 

If in addition / £ M[x], the sequence z = (z a ) is obtained as a linear combination 
of (y Q ). Indeed if /(x) = ^ /^x" then 

z Q = X] f0y a +0> VaeN", 
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and so in (|3.4[) . Md(z) is nothing less than the localizing matrix M<j(/y) associated 
with y = (y a ) and / £ K[x], defined in ([231), and (0IU) reads M d (/y) b for all 
d = 0,l,... 

Finally, if / is concave then / > on K implies / > on co (K), and so the only 
i/part of (b) also follows. □ 

Therefore, to check whether a polynomial / E R[x] is nonnegative on K, it 
suffices to check if every element of the countable family of real symmetric matrices 
(Md(/y)), d E N, is positive semidefinite. 



Remark 3.3. An informal alternative proof of Theorem 13.21 which does not use 
Lemma 13.11 is as follows. If / is not nonnegative on K there exists a E K such 
that /(a) < 0, and so as K is compact, there is a continuous function, e.g, x 
h(x) := exp(— c||x — a|| 2 ) close to 1 in some open neighborhood B(a, S) of a, and 
very small in the rest of K. By the Stone- Weierstrass's theorem, one may choose 
h to be a polynomial. Next, the complement B(a, S) c (= R" \ B(a, 5)) of B(a, 5) is 
closed, and so K n B(a, <5) c is a closed set contained in K (hence smaller than K). 
Therefore /i(B(a, 6)) > because otherwise /i(K n B(a, i5) c ) = /i(K) which would 
imply that K n B(a, S) c is a support of fi smaller than K, in contradiction with 
supp/i = K. Hence, we would get the contradiction 



Jh 2 fdfi^ h(a) 2 f(a)fx(B( a ,6)) < 0. 



However, in the non compact case described in the next section, this argument is 
not valid. 

3.2. The non-compact case. We now consider the more delicate case where K 
is a closed set of M n , not necessarily compact. To handle arbitrary non compact 
sets K and arbitrary continuous functions /, we need a reference measure p, with 
supp/i = K and with nice properties so that integrals such as J g 2 fd[i, g E M[x], 
are well-behaved. 

So, let ip be an arbitrary finite Borel measure on R n whose support is exactly 
K, and let (x be the finite Borel measure defined by: 



(3.7) 



H{B) := J exp ^-J^mI dtp{x), VBe6(K"). 



Observe that supp/i = K and /x satisfies Carleman's condition (|2.2[) . Indeed, let 
z = (z a ), a E N™, be the sequence of moments of /i. Then for every i — 1, . . . , n, 
and every k = 0, 1, . . ., using x\ k < (2fe)! exp \xi\, 

(3.8) L z (x 2k )= [ xf dfi(x) < (2k)\ ( e^d/^x) < (2fc)!p(K) =:{2k)\M. 
Therefore for every i = 1, . . . , n, using (2fc)! < (2k) 2k for every k, yields 



co co co M -i/2k 

Y^L^xf)- 1 ^ > VA/- 1 / 2fc ((2fc)!)- 1 ^> V^-— = +oo, 

fc=i fc=i fc=i 

i.e., (|2.2[) holds. Notice also that all the moments of fi (defined in (|3.7jl ) are finite, 
and so every polynomial is /i-integrable. 
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Theorem 3.4. Let K C R n be closed and let tp be an arbitrary finite Borel measure 
whose support is exactly K. Let f be a continuous function on R n . If f € R[x] 
(i.e., f is a polynomial) let /i be as in \3. 7j ) whereas if f is not a polynomial let /i 
be defined by 

(3.9) v(B) := / ^Pi-SJ^^(x), VBee(l"). 

TTien (a,) and (b) of Theorem \3.2\ hold. 
For a detailed proof see 

It is important to emphasize that in Theorem 13.21 and 13.41 the set K is an 
arbitrary closed set of K" , and to the best of our knowledge, the characterization of 
nonnegativity of / in terms of positive definiteness of the moment matrices Mj(z) 
is new. But of course, this characterization becomes even more interesting when 
one knows how the compute the moment sequence z = (z a ), a € N n , which is 
possible in a few special cases only. 

Important particular cases of nice such sets K include boxes, hyper rectangles, 
ellipsoids, and simplices in the compact case, and the positive orthant, or the whole 
space R n in the non compact case. For instance, for the whole space K = R n one 
may choose for /j, in (|3.7[) the multivariate Gaussian (or normal) probability measure 

fi(B) := (2tt)-"/ 2 exp(-i||x|| 2 ) dx, B e B(R n ), 

which the rt-times product of the one-dimensional normal distribution 

Hi{B) := -L / exp(-a;?/2) <fci, B £ B(R), 
v27r Jb 

whose moments are all easily available in closed form. In Theorem 13.41 this corre- 
sponds to the choice 

(3.10) <p{B) = (2^)-/ 2 / ex P(HN 2 / 2 ) dX; B G B{R ny 

Jb cxp(-^ l=1 \Xi\) 

When K is the positive orthant E™ one may choose for /j, the exponential probability 
measure 

(3.11) p{B) := / cxp(- Y\xi)dx, B e B(M"), 

Jb i=i 

which the n-times product of the one-dimensional exponential distribution 

m{B) := / exp(~x,)dxi, B € B(R+), 
Jb 



whose moments are also easily available in closed form. In Theorem 13.41 this corre- 
sponds to the choice 

/n 
exp(- J^xjdx, B e B(M"). 
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3.3. The cone of nonnegative polynomials. The convex cone Cd C R[x]2d of 
nonnegative polynomials of degree at most 2d (a nonnegative polynomial has nec- 
essarily even degree) is much harder to characterize than its subcone E[x]<j of sums 
of squares. Indeed, while the latter has a simple semidefinite representation with 
lifting (i.e. E[x]<j is the projection in R s ( 2d ' of a spectrahedroijj in a higher di- 
mensional space), so far there is no such simple representation for the former. In 
addition, when d is fixed, Blekherman [4] has shown that after proper normaliza- 
tion, the "gap" between Cd and S[x]d increases unboundedly with the number of 
variables. 

We next provide a convergent hierarchy of (outer) semidefinite approximations 
(C^), k £ N, of Cd where each C\ has a semidefinite representation with no lift- 
ing (i.e., no projection is needed and is a spectrahedron). To the best of our 
knowledge, this is the first result of this kind. 

Recall that with every / £ R[x]d is associated its vector of coefficients f = (f a ), 
a £ Nj, in the canonical basis of monomials, and conversely, with every f £ R s ( d ) 
is associated a polynomial / £ M.[x]d with vector of coefficients f = (f a ) in the 
canonical basis. Recall that for every k = 1, . . ., 

7 =J- r X r e -* 2 /*d X - I k ° lf ^ = 2fc + 1 ' 



as 



12k = (2fc - l)72(fc-i) for every k > 1. 



Corollary 3.5. Let fj, be the probability measure on R ra which is the n-times product 
of the normal distribution on R, and so with moments y = (y a ), ct £ N n , 



(3.12) y. 



a = f x a dn = f] f x a *e- x2/2 dx) , Vae 

JR" i=1 V V27T J-oo / 



For every k £ N, let C k d := {f £ W ( - 2d '> : M k (fy) y 0}, where M k (fy) is the 
localizing matrix in V2.3\) associated with y and f £ R[x] 2( j. Each C\ is a closed 
convex cone and a spectrahedron. 

Then D C\ ■ ■ ■ D C\ ■ ■ ■ D Cd and f £ Cd if and only if its vector of coefficients 
f £ R s ( 2d ) satisfies feCj, for every k — 0,1, 

Proof. Following its definition (|2.3[) . all entries of the localizing matrix Mfe(/y) are 
linear in f £ R s ( 2d ), and so M fe (/y) h is an LMI. Therefore C^ is a spectrahedron 
and a closed convex cone. Next, let K := R" and let /i be as in Corollary 13.51 and so 
of the form (|3.7p with (p as in (|3.10[) . Then /i satisfies Carleman's condition (|2.2p . 
Hence, by Theorem 13.41 with K = R", / is nonnegative on K if and only if p.4[) 
holds, which is equivalent to stating that M^(/y) ^0, k = 0, 1, . . ., which in turn 
is equivalent to stating that f £ C^, k — 0,1, . . . □ 

So the nested sequence of convex cones C° D C\ ■ ■ ■ D Cd defines arbitrary close 
outer approximations of Cd- In fact C\^ =l fi\ is closed and Cd = n^ C^. It is worth 
emphasizing that each C^ is a spectrahedron with no lifting, that is, C% is defined 
solely in terms of the vector of coefficients f with no additional variable (i.e., no 
projection is needed). 



A spectrahedron is the intersection of the cone of positive semidefinite matrices with an 
affinc-lincar space. Its algebraic representation is called a Linear Matrix Inequality (LMI). 
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For instance, the first approximation C° is just the set {f G K s ( 2d ) ; J fd/i > 0}, 
which is a half-space of R s ( 2d ). And with n — 2, 



f e 



fdfi J x 1 fdfi J x 2 fdfi 
xifdn / xlfdfi / x 1 x 2 fd/j, 



x 2 fdfi \ x\x 2 fd\i \ x 2 fd/i 



> 



or, equivalently, C\ is the convex basic semi-algebraic set: 
jf G K s(2d) : y fdfi > 

*?/d^ (y /d/i) > (y *</d/i) , «= 1,2 
y %ifdf^j ^y xifdiA > ^y x 1 x 2 fd^ 

fd[i\ xffdfij ^ xlfdfxj - (Jxix 2 fdi^j 
J xifdfij (Jxlfdij\ - ^ x 2 fdf?j (Jxlfdfi 



xxfd/j, 



x 2 fdfi 



X\x 2 fd[i I > 



where we have just expressed the nonnegativity of all principal minors of Mi(/y). 

A very similar result holds for the convex cone Cd(K) of polynomials of degree 
at most d, nonnegative on a closed set K C K™. 



Corollary 3.6. Let K C K" be a closed set and let /i be defined in \3. 7| ) where ip 
is a an arbitrary finite Borel measure whose support is exactly K. 

For every k G N, letC^(K) := {f G : M k (fy) t 0}, where M fc (/y) is i/ie 

localizing matrix in [2.3]) associated with y and f G M[x]rf. iSac/i C^'(K) is a closed 
convex cone and a spectrahedron. 

Then C°(K) D C](K) • • O Cj(K) • • O C d (K) and / € C d {K) if and only if its 
vector of coefficients f G R s ( d ) satisfies f € C^(K), /or every k = 0, 1, . . .. 



The proof which mimicks that of Corollarv l3.5l is omitted. Of course, for practical 
computation, one is restricted to sets K where one may compute effectively the 
moments of the measure /i. An example of such a set K is the positive orthant, 
in which case one may choose the measure [i in (|3.1ip for which all moments are 
explicitly available. For compact sets K let us mention balls, boxes, ellipsoids, 
and simplices. But again, any compact set where one knows how to compute all 
moments of some measure with support exactly K, is fine. 

To the best of our knowledge this is the first characterization of an outer ap- 
proximation of the cone C^(K) in a relatively general context. Indeed, for the basic 
semi-algebraic set 



(3.13) 



K 



{xG 



9j{x) > 0, j = l,...,m}, 
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Stengle's Nichtnegativstellensatz [22] states that / G K[x] is nonnegative on K if 
and only if 

(3.14) Pf = f s +q, 

for some integer s and polynomials p,q G P(g), where P(g) is the preorderin^ 
associated with the g/s. In addition, there exist bounds on the integer s and the 
degree of the s.o.s. weights in the definition of p, q G P(g), so that in principle, 
when / is known, checking whether / > on K reduces to solving a single SDP 
to compute h,p in the nonnegativity certificate (|3.14|) . However, the size of this 
SDP is potentially huge and makes it unpractical. Moreover, the representation of 
Cd(K) in (|3.14[) is not convex in the vector of coefficients of / because it involves 
f 2s as well as the product pf. 

Remark 3.7. If in Corollary 13.61 one replaces the finite-dimensional convex cone 
Crf(K) C R[x]d with the infinite-dimensional convex cone C(K) C K[x] of all poly- 
nomials nonnegative on K, and C*(K) C M[x] d with C k (K) = {f G R s ( 2fe ) : 
Mfc(/y) y 0}, then the nested sequence of (increasing but finite-dimensional) 
convex cones C k (K), k G N, provides finite-dimensional approximations of C(K). 

4. Application to polynomial optimization 
Consider the polynomial optimization problem 

(4.1) P: f* = inf{/(x) : xGK}, 

X 

where K C 1" is closed and / G R[x]. 

If K is compact let /i be a finite Borel measure with supp fi = K and if K is not 
compact, let ip be an arbitrary finite Borel measure with supp (p — K and let /i be as 
in (|3.7p . In both cases, the sequence of moments y = (y a ), a G N n , is well-defined, 
and we assume that y a is available or can be computed, for every a G N™. 

Consider the sequence of semidefinite programs: 

(4.2) A d = sup{A : M d (/ A y)^0} 

age 

where fx G R[x] is the polynomial x M> /(x) — A. Notice that (|4.2I) has only one 
variable! 

Theorem 4.1. Consider the hierarchy of semidefinite programs indexed by 

d G N. Then: 

(i) {4-*fy has an optimal solution Xd > /* for every d G N™. 

(ii) The sequence (\d), d G N, is monotone nonincreasing and Xd i f* as d — > oo. 

Proof, (i) Since / — /* > on K, by Theorem l3.21 A := /* is a feasible solution of 
(|4.2j) for every d. Hence Xd > f* for every deN. Next, let d G N be fixed, and let 
A be an arbitrary feasible solution of (j4.2|) . From the condition M^/xy) h 0, the 
diagonal entry M.d(f\ y)(l, 1) must be nonnegative, i.e., Ayo < J2 a faVa, and so, 
as we maximize and yo > 0, (|4. 2[) must have an optimal solution Ad. 

(ii) Obviously Ad < A m whenever d > m, because Md{f\y) h implies 
M m (/Ay) >z 0. Therefore, the sequence (Ad), d G N, is monotone nonincreas- 
ing and being bounded below by /*, converges to A* > /*. Next, suppose that 



The preordering P{g) associated with the gj's is the set of polynomials of the form 
E o6 {o,l) m ° a Bi 1 ■■ ■ 9m m where a a e E[x] for each a. 
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A* > /*; fix k e N, arbitrary. The convergence Ad i A* implies M.k(f\* y) h 0. 
As k was arbitrary, we obtain that ~M.d(f\* y) h for every d E N. But then by 
Theorem 13.21 or Theorem 13. 4[ / — A* > on K, and so A* < /*, in contradiction 
with A* > /*. Therefore A* = /*. □ 



For each d E N, the semidefinitc program (|4.2[) provides an upper bound on the 
optimal value /* only. We next show that the dual contains some information on 
global minimizers, at least when d is sufficiently large. 

4.1. Duality. Let S d be the space of real symmetric s(d) x s(d) matrices. One 
may write the semidefinitc program (|4.2I) as 

(4.3) X d = sup{A : AM d ( y ) * M d (fy)}, 

x 

which in fact is a generalized eigenvalue problem for the pair of matrices M d (y) 
and M d (/y). Its dual is the semidefinite program 



inf { (X,M d (/y)) : (X,M d (y)) = 1; X h 0}, 



or, equivalently, 



(4.4) A d = inf jy fadfi: J adfj. = 1; a G S[x] d i 



So the dual problem (|4.4j) is to find a sum of squares polynomial a of degree at 
most 2d (normalized to satisfy J adfi = 1) that minimizes the integral J fcrdfi, and 
a simple interpretation of (|4.4I) is as follows: 

With M(K) being the space of Borel probability measures on K, we know that 

f* = inf / fdip. Next, let Md(ff) C M(K) be the space of probability mea- 
v eAf(K) J K 

sures on K which have a density a € S[x] d with respect to \x. Then (I4.4p reads 



inf / fdip, which clearly shows why one obtains an upper bound on /*. In- 

deed, instead of searching in M (K) one searches in its subset M d (/i). What is not 
obvious at all is whether the obtained upper bound obtained in (|4.4|) converges to 
/* when the degree of a E S[x] d is allowed to increase! 

Theorem 4.2. Suppose that f* > —oo and K has nonempty interior. Then : 

(a) There is no duality gap between {4-2ty and |^.^[ ) and J^.^[ ) has an optimal 

solution a* E S[x](j which satisfies / (/(x) — A d ) <r*(x)<i/i(x) = 0. 

Jk 

(b) // K is convex and f is convex, let s.* d := J xer*(x)<i^i(x). Then x d 6 K 
and /* < /(x d ) < A d , so i/iat /(xj) — > /* as d — > oo. Moreover, if the set 
{x G K : /(x) < /o} is compact for some fo > f* , then any accumulation point 
x* E K of f/ie sequence (x d ), d € N, is a minimizer of problem that is, 
/(x*) = /*. 

Proof, (a) Any scalar A < /* is a feasible solution of (|4.2p and in addition, M d ((/ — 
A) y) ^0 because since K has nonempty interior and / — A > on K, 

(g,M d ((/-A)y)g) = / (/(x) - A) ff (x)V(dx) > 0, V 9 eK[x] d . 
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But this means that Slater's condition^ holds for (|4.2p which in turn implies that 
there is no duality gap and (|4.4[) has an optimal solution a* € £[x]d; see e.g. [23] . 
And so, 

/ (/(x)-A d )<7*(x)d/i(x) = / fa*d(i-\ d = 0. 

(b) Let v be the Borel probability measure on K defined by v(B) = J B a*dfi, B € B. 
As / is convex, by Jensen's inequality (see e.g. McShane [15]). 

jj°*dn = Jjdv> fUxdA = fiO- 

In addition, if K is convex then 6 K and so, A^ > /(xj) > /*. Finally if for 
some /o > /*, the set % := {x € K : /(x) < /o} is compact, and since A^ — > /*, 
then f(xj) < fo for <i sufficiently large, i.e., x. d e "H for sufficiently large d By 
compactness there is a subsequence (dg), I £ N, and a point x* 6 K such that 
x^ — > x* as £ — > oo. Continuity of / combined with the convergence /(x^) — > /* 
yields /(xj ) — > /(x*) = /* as £ — >• oo. As the convergent subsequence (x5 ) was 
arbitrary, the proof is complete. □ 

So in case where / is a convex polynomial and K is a convex set, Theorem 14.21 
provides a means of approximating not only the optimal value /* , but also a global 
minimizcr x* € K. 

In the more subtle nonconvex case, one still can obtain some information on 
global minimizers from an optimal solution a* E S[x]d of (|4.4p . Let e > be fixed, 
and suppose that d is large enough so that f*<Xd<f*+ £• Then, by Theorem 
IHa), 

/ (/(x)-/>*(x)d/i(x) =X d -f* <e. 

As / — /* > on K, necessarily the measure dv = a* dp, gives very small weight to 
regions of K where /(x) is significantly larger than /*. For instance, if e = 10~ 2 
and A := {x e K : /(x) > /* + 1}, then u(A) < 10~ 2 , i.e., the set A contributes to 
less than 1% of the total mass of v. So if /x is uniformly distributed on K (which is a 
reasonable choice if one has to compute all moments of fi) then a simple inspection 
of the values of <r*(x) provides some rough indication on where (in K) /(x) is close 
to /*. 

The interpretation (|4.4j) of the dual shows that in general the monotone conver- 
gence is only asymptotic and cannot be finite. Indeed if K has a nonempty interior 
then the probability measure dv — ad\x cannot be a Dirac measure at any global 
minimizer x* £ K. An exception is the discrete case, i.e., when K is a finite number 
of points, like in e.g. 0/1 programs. Indeed we get: 

Corollary 4.3. Let KcB" be a discrete set (x(fc)) C R", k € J, and let \i be the 
'probability measure uniformly distributed in K, i.e., 

1 s 

m = -z2 6 x(fc), 
fc=i 

where s — \ J\ and i5 x denote the Dirac measure at the point x. Then the optimal 
value Xd of \4-2^ satisfies \ c i = f* for some integer d. 

^For an optimization problem inf x {/o(x) : /j(x) > 0, j = 1, . . . ,m}, Slater's condition states 
that there exists xq such that /j(xo) > for every j = 1, ,m. 
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Proof. Let x* = x(j*) (for some j* € J) be the global minimizer of min{/(x) : 
x € K}. For each k = 1, . . . , s there exists a polynomial gfc € R[x] such that 
<7/c(x(j)) = 8k=j for every j = 1, . . . , s (where 5k=j denotes the Kronecker symbol). 
The polynomials (%) are called interpolation polynomials. So let a* := sg|, € E[x], 
so that 

/ /(x)<7*(x)d/x(x) = /(x(f )) = /* and f a*dfi = 1. 

Hence as soon as d > deg <?.,■», <r* € E[x]d is a feasible solution of (14.41) . and so 
from /* = j fa* dpi > X* d > \d > f* we deduce that = = /*, the desired 
result. □ 

There are several interesting cases where the above described methodology can 
apply, i.e., cases where y can be obtained either explicitly in closed form or numer- 
ically. In particular, when K is cither: 

• A box B := nr=i[ a »'^] c wrt h A 4 being the normalized Lebesgue 
measure on B. The sequence y = (y a ) is trivial to obtain in closed form. 

• The discrete set { — 1, 1}™ with /i being uniformly distributed and normal- 
ized. Again the sequence y = (y a ) is trivial to obtain in closed form. Notice 
that in this case we obtain a new hierarchy of semidefinite relaxations (with 
only one variable) for the celebrated MAXCUT problem (and any nonlinear 
0/1 program). 

• The unit Euclidean ball B:={x : ||x|| 2 < 1} with /i uniformly distributed, 
and similarly the unit sphere S := {x : ||x|| 2 = 1}, with \i being the 
rotation invariant probability measure on S. In both cases the moments 
y = (ya) are obtained easily. 

• A simplex A c K™, in which case if one takes \x as the Lebesgue measure 
then all moments of fi can be computed numerically. In particular, with d 
fixed, this computation can be done in time polynomial time. See e.g. the 
recent work of [2]. 

• The whole space M™ in which case /i may be chosen to be the product 
measure <8>™ =1 z^ with each being the normal distribution. Observe that 
one then obtains a new hierarchy of semidefinite approximations (upper 
bounds) for unconstrained global optimization. The corresponding mono- 
tone sequence of upper bounds converges to /* no matter if the problem 
has a global minimizer or not. This may be an alternative and/or a com- 
plement to the recent convex relaxations provided in Schweighofer [21] and 
Ha and Vui [8] which also work when /* is not attained, and provide a 
convergent sequence of lower bounds. 

• The positive orthant R" , in which case /i may be chosen to be the product 
measure ^—l^i with each Vi being the exponential distribution Vi{B) — 
f M+nB e~ x dx, B 6 B. In particular if x h-> f(x) := x T Ax where A e S n , 
then one obtains a hierarchy of numerical tests to check whether A is a 
copositive matrix. Indeed, if is an optimal solution of (|4.3|) then A is 
copositive if and only if Ad > for all d € N. Notice that we also obtain a 
hierarchy of outer approximations (COP^) C S n of the cone COP ofnxn 
copositive matrices. Indeed, for every Ag5„, let /a be the quadratic form 
x i — ^ /a(x) := x T Ax. Then, for every d, the set 

COP d := {AeS n : M d (/ A y) ^0} 
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is a convex cone denned only in terms of the coefficients of the matrix A. 
ft is even a spectrahedron since M^(/a y) is a linear matrix inequality in 
the coefficients of A. And in view of Theorem l3.2f a). COP = C\ deN COP^- 

4.2. Examples. In this section we provide three simple examples to illustrate the 
above methodology. 

Example 1. Consider the global minimization on K = of the Motzkin-like 
polynomial x H> /(x) = x\x\ {x\ + a; 2 — 1) whose global minimum is /* = —1/27 s» 
—0.037, attained at (sc*,^) = (±-\/f /3, ±^/l/3). Choose for fi the probability 
measure n(B) := J B e~ xl ~ X2 dx, B £ B(M.+), for which the sequence of moments 
y = (yij), i,j £ N, is easy to obtain. Namely jjij = for every i,j £ N. Then 
the semidchnite relaxations (|4.2p yield Ao = 92, Ai = 1.5097, and A14 = —0.0113, 
showing a significant and rapid decrease in first iterations with a long tail close 
to /*, illustrated in Figure [TJ Then after d = 14, one encounters some numerical 
problems and we cannot trust the results anymore. 

If we now minimize the same polynomial / on the box [0, l] 2 , one choose for /i 
the probability uniformly distributed on [0, l] 2 , whose moments y = i,j £ N, 

are also easily obtained by = (i + 1) _1 0' + l) -1 . Then one obtains Ao = 0.222, 
Ai = —0.055, and A10 = —0.0311, showing again a rapid decrease in first iterations 
with a long tail close to /*, illustrated in Figure [2j 




Figure 1 . Minimizing the Motzkin-like polynomial in K 2 



Example 2. Still on K = ffi^_ , consider the global minimization of the polynomial 
x i— > x\ + (1 — X1X2) 2 whose global minimum /* = is not attained. Again, 
choose for \i the probability measure (J,(B) := f B c~ Xl ~ X2 dx, B £ ^(K. 2 ^). Then the 
semidefinite relaxations (|4.2p yield Ao = 5, Ai = 1.9187 and A15 = 0.4795, showing 
again a significant and rapid decrease in first iterations with a long tail close to /*, 
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0.03 




Figure 2. Minimizing the Motzkin-like polynomial in [0, l] 2 

illustrated in Figure [3] numerical problems occur after d — 15. However, this kind 
of problems where the global minimum /* is not attained, is notoriously difficult. 
Even the semidefinite relaxations defined in [8] (which provide lower bounds on /*) 
and especially devised for such problems, encounter numerical difficulties; see [8j 
Example 4.8]. 




i 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 



FIGURE 3. Minimizing x\ + (1 - xix 2 ) 2 on R 2 
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Example 3. The following example illustrates the duality results of Section §4.11 
The univariate polynomial x ^ f(x) := 0.375 — 5x + 21x 2 — 32x 3 + 16a; displayed in 
FigHhas two global minima at x\ = 0.1939 and x* 2 = 0.8062, with /* = -0.0156. 
In Fig [5] is plotted the sequence of upper bounds Ad — > /* as d — > oo, with again 
a rapid decrease in first iterations. One has plotted in Fig [B] the s.o.s. polynomial 
x i y <i(x), optimal solution of (|4.4p with d = 10, associated with the probability 
density a(x)dx as explained in jj4.il As expected, two peaks appear at the points 

Xi ~ X^ , z 1, 2. 



0.4 




Figure 4. f(x) = 0.375 - bx + 21x 2 - 32x a + 16x 4 on [0, 1] 

Example 4. We finally consider a discrete optimization problem, namely the cel- 
ebrated MAXCUT problem 

/* = min{x T Qx : x€ {-1,1}"}, 

X 

where Q = (Qij) £ M. nxn is a real symmetric matrix whose all diagonal elements 
vanish. The measure /i is uniformly distributed on { — 1, 1}™ so that its moments 
are readily available. We first consider the equal weights case, i.e., Qij = 1/2 for 
all with i ^ j in which case /* = — [n/2\. With n = 11 the successive values 
for Ad, d < 4, are displayed in Tableland A4 is relatively close to /*. Next we 
have generated five random instances of MAXCUT with n = 11 but Qij = with 
probability 1/2, and if Qij ^ it is randomly generated using the Matlab "rand" 
function. The successive values of Ad, d < 4, are displayed in Table [5J and again, 
A 4 is quite close to /"Q. 

The above examples seem to indicate that even though one chooses a measure /i 
uniformly distributed on K, one obtains a rapid decrease in the first iterations and 

^The optimal value /* has been computed using the GloptiPoly software I10| dedicated to 
solving the Generalized Problem of Moments 
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0.08 




FIGURE 5. Minimizing 0.375 - 5a; + 21x 2 - 32a; 3 + 16a; 4 on [0, 1] 




FIGURE 6. The probability density a{x)dx on [0, 1] 

then a slow convergence close to /*. If on the one hand the convergence to /* is 
likely to be slow, on the other hand, one has to solve semidefmite programs (|4.2[) 
with only one variable! In fact solving the semidefinite program (|4.3j) is computing 
the smallest generalized eigenvalue associated with the pair of real symmetric ma- 
trices (Md(/y),Md(y)), for which specialized codes are available (instead of using 
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d 


d = 


d=l 


d = 2 


d = 3 


d = A 


r 







-1 


-2.662 -3.22 


-4 


-5 



Table 1. MAXCUT: n = 11; Q(i, j) = 1 for all i ^ j. 



d 


Ao 


Ai 


A 2 


A 3 


A 4 


/* 


Exl 





-1.928 


-3.748 


-5.22 


-6.37 


-7.946 


Ex2 





-1.56 


-3.103 


-4.314 


-5.282 


-6.863 


Ex3 





-1.910 


-3.694 


-5.078 


-6.161 


-8.032 


Ex4 





-2.164 


-4.1664 


-5.7971 


-7.06 


-9.198 


Ex5 





-1.825 


-3.560 


-4.945 


-5.924 


-7.467 



Table 2. MAXCUT: n = 11; Q random. 



a solver for semidefinite programs). However, one has to remember that the choice 
is limited to measures fi with supp \x — K and whose moments are available or easy 
to compute. Hence, the present methodology is so far limited to simple sets K as 
described before. Finally, analyzing how the convergence to /* depends on /i is 
beyond the scope of the present paper and is a topic of further research. 

4.3. Discusssion. In nonlinear programming, sequences of upper bounds on the 
global minimum /* are usually obtained from feasible points x £ K, e.g., via some 
(local) minimization algorithm. But for non convex problems, providing a sequence 
of upper bounds that converges to the global minimum /* is in general impossible 
unless one computes points on a grid whose mesh size tends to zero. In the above 
methodology one provides a monotone nonincreasing sequence of upper bounds 
converging to /* for polynomial optimization problems on sets K, non necessarily 
compact but such that one may compute all moments of some finite Borel measure 
fi with supp \x = K. In fact, if there are only finitely many (say up to order 2d) 
moments available then one obtains a finite sequence of upper bounds. 

In contrast to the hierarchy of semidefinite relaxations in e.g. [12l US] which 
provide lower bounds converging to /* when K is a compact basic semi- algebraic 
set, the convergence of the upper bounds to /* is only asymptotic and never finite, 
except when K is a discrete set. However, and even if we expect the convergence 
to be rather slow when close to /* , to our knowledge it is the first approach of this 
kind, and in a few iterations one may obtain upper bounds which (even if crude) 
complements the lower bounds obtained in [12] (in the compact case). 

Also note that to solve (|4.3|) several improvements are possible. For instance, we 
have already mentioned that it could be solved via specialized packages for general- 
ized eigenvalue problems. Next, if instead of using the canonical basis of monomial 
(x Q ), one now expresses the moment matrix Mrf(y) with rows and columns indexed 
in the basis of polynomials (p a ) C R[x] (up to degree d) orthogonal with respect 
to n, then M^y) becomes the identity matrix. And so problem (|4.3j) reduces to 
a standard eigenvalue problem, namely that of computing the smallest eigenvalue 
of the (real and symmetric) localizing matrix Md(/y) (expressed in the basis of 



family of univariate polynomials (pj.) C M[x] is orthogonal with respect to a finite measure 
fi on M if f piPkdfi = For extensions to the multivariate case see e.g. [6] [9]. 
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orthogonal polynomials)! And it turns out that computing the orthogonal poly- 
nomials is easy once the moment matrix M<j(y) is available, since they can be 
obtained via computing certain determinants, as explained in e.g. (5J[9]. 

Inverse problem from moments. Finally, observe that the above methodology 
perfectly fits inverse problems from moments, where precisely some Borel measure 
fj, is known only from its moments (via some measurement device), and one wishes 
to recover (or approximately recover) its support K from the known moments; see 
e.g. the work of Cuyt et al. [5] and the many references therein. Hence if / € M[x] 
is fixed then by definition / — /* > (on K) provides a strong valid (polynomial) 
inequality for the unknown set K. So computing an optimal solution Xd of (|4.2[) for 
d sufficiently large, will provide an almost- valid polynomial inequality / — A<j > 
for K. One may even let / G M[x]<2 be unknown and search for the "best" valid 
inequality /(x) — /* > where / varies in some family (e.g. linear or quadratic 
polynomials) and minimizes ome appropriate (linear or convex) objective function 
of its vector of coefficcnts f . 

5. Conclusion 

In this paper we have presented a new characterization of nonnegativity on a 
closed set K which is based on the knowledge of a single finite Borel measure \x 
with supp (i = K. It permits to obtain a hierarchy of spectrahedra which provides 
a nested sequence of outer approximations of the convex cone of polynomials of 
degree at most d, nonnegative on K. When used in polynomial optimization for 
certain "simple sets" K, one obtains a hierarchy of semidefinite approximations 
(with only one variable) which provides a nonincreasing sequence of upper bounds 
converging to the global optimum, hence a complement to the sequence of upper 
bounds provided by the hierarchy of semidefinite relaxations defined in e.g. p~2l [T3] 
when K is compact and basic semi- algebraic. A topic of further investigation is to 
analyze the efficiency of such an approach on a sample of optimization problems 
on simple closed sets like the whole space W 1 , the positive orthant R™ , a box, a 
simplex, or an ellipsoid, as well as for some inverse problems from moments. 

6. APPENDIX 

Proof of Theorem 13.41 

Proof. The only z/part is exactly the same as in the proof of Theorem l3.2l For the 
i/part, let z = (z a ) be the sequence defined in (|3.5|) . The sequence z is well defined 
because x x Q /(x) is /i-integrable for all a £ N n . Indeed, if / is a polynomial (so 
that \i is defined in (|3.7|) ) we have seen that all moments of fi are finite and since 
x Q /(x) is a polynomial the result follows. If / is not a polynomial (so that /i is 
defined in ^9\) then 

f |x a /(x)|d/x(x) < f |x Q | exp(- V \ Xi \) <Mx) <<p(M n ) 17^!, 

where we have used that Ix"*] < c^Iexp |a^|, and |/|/(1 + f 2 ) < 1 for all x. As in 
the proof of Theorem l3.2[ the set function B M> v(B) := J B fd/i, B G B, is a signed 
Borel measure because again v can be written as the difference v\ — V2 of the two 
positive Borel measures z^i,^2 in (|3.1[) - (|3.2[) . With same majorizations as above, 
both v\ and Vi are finite Borel measures and so v is a finite signed Borel measure. 
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The same arguments as in the proof of Theorem 13.21 show that Mrf(z) >z for 
every d £ N. Next, the sequence z satisfies the generalized Carleman's condition 



Indeed, first consider the case where / is a polynomial (and so /i is as in (13. 71) ). 
Let 1 < i < n be fixed arbitrary, and let 2s > deg/. Observe that whenever \a\ < k, 
|x| Q < |xj| fe on the subset Wj := {x £ E" \ [-1, 1]" : \xj\ = max, \ Xi \}. And so, 
|/(x)| < H/llx | Xj | 2s for all x e W 3 (and where Wfh := £ Q !/«!)■ Hence, 



L z (xf) = / /(x)^ fe ^(x) 
Jk 



< / |/(x)|xfdMx) + ||/||iE / ^f +S) ^(x) 

JKn[-i,i]» J=1 JKnff,, 

(6.1) < H/ll !/i(K) + Mn||/|| ! (2(fc + a))! < 2Mn||/||i (2(As + «))!, 

where Af is as in (|3.8[) and assuming with no loss of generality that /i(K) < 
Mn(2(k + s))l (otherwise rescale </?). And so we have 



Lz(xfT 1/2fe > (2Mn\\f\\ 1 )- 1 / 2k ((2(fc + S ))!)- 1 / 2 ( fc+s )) 
> 1 ((2(fc + S ))!)- 1 / 2 ( fe + s )) 



0+s)/fc 



(fe+s)/fe 



> - 



I / 1 \ (k+s)/k 



2 \2(k + s) 

where k > ko is sufficiently large so that {2Mn\\f\\i)~ 1 / 2k > 1/2. Therefore, 

oo oo , s (fe+s)/fc 

fe=l fc=fc Q v v ;/ 

where the last equality follows from X^fcLi^^) 1 = +oo. Indeed, ( 2 (k+s) Y k+S ^ k = 
( ■2(k+s) )( 2(k+ s ) ) s/k and (. 2(k+ s ) y /k - 1 / 2 whenever k is sufficiently large, say k > 
k\. Hence the sequence z satisfies Carleman's condition (|2.2[) . 
If / is not a polynomial then fi is as in (|3.9j) and so 

L,{xf) = J,f*|kM/(x)#(x) 

„ 2 fc ex p(-Er=i N) 



(6.2) 



< y *f \Xp - 1/Wldy(x) 

< y x 2fc exp ("XJ NJ ^« ^ ( 2fc ) !M > 



where we have used that |/|/(1 + / 2 ) < 1 for all x, and x 2k < (2fc)!exp \x{\. And 
so again, the sequence z satisfies Carleman's condition (I2.2j) . 

Next, as Md(z) >; for every d £ N, by Proposition l2.11 z is the moment sequence 
of a measure ip on R" and -0 is determinate. In addition, from the definition p.5[) 
of i/ and z, we have 

(6.3) f x Q ^(x)=z Q = f x Q d^(x), VaeN" 
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But then using v = v\ — V2 in (|3.1[) - (|3.2[) . (|6.3|) reads 



(6.4) f x a d(^ + i/ 2 )(x) = f x a di*(x), 



Va G N". 




Let v = (v a ), a G N™, be the sequence of moments associated with v\. Of course, 
Md(v) >: for all del Next, 



and so, depending on whether / is a polynomial or not, we obtain L v (xf k ) < 
2Mn||/||i(2(fe + s))! as in dHUJ) or L v (xf k ) < (2k)\M as in In both cases the 

sequence v satisfies Carleman's condition (|2.2[) and since Md(v) >: for all d G N, 
by Proposition ^. 1[ i/i is moment determinate. But then (|6.4I) yields = ip + 1/2, 
or equivalently, V- 1 = *T — ^2 (= that is, v is a positive measure. The rest of the 
proof is exactly the same as for proof of Theorem 13.21 □ 
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